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We examine the 1/N expansion, where N is the number of two-component Dirac fermions, for 
Coulomb interactions in graphene with a gap of magnitude A = 2m. We find that for Net S> 1, 
where a is graphene's "fine structure constant", there is a crossover as a function of distance r 
from the usual 3D Coulomb law, V(r) ~ 1/r, to a 2D Coulomb interaction, V(r) ~ \n(Na/mr), 
for m _1 <r< m _1 Na/6. This effect reflects the weak "confinement" of the electric field in the 
graphene plane. The crossover also leads to unusual renormalization of the quasiparticle velocity 
and gap at low momenta. We also discuss the differences between the interaction potential in 
gapped graphene and usual QED for different coupling regimes. 
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I. INTRODUCTION 

The physics of graphene, a two-dimensional (2D) al- 
lotrope of carbon, presents a unique opportunity to ex- 
plore properties of gapless (massless) Dirac fermions in 
a solid-state context^. Because graphene is a truly 2D 
material from electronic point of view, it is also an inter- 
esting system where one can explore the important issue 
of electron-electron interactions. The interaction prob- 
lem is a central one for the physics of low-dimensional 
electron systems. Although the electron-electron inter- 
actions in graphene are not internally screened, that is, 
they remain long-range and decay like 1/r (the electric 
field lines propagate in 3D away from the graphene plane 
and one has the ordinary Coulomb law), there is little 
experimental evidence, if any, that electron-electron in- 
teractions play a role in graphene physics. It is possi- 
ble that the interactions between the fermions are actu- 
ally dielectrically screened by the presence of substrates, 
onto which graphene is deposited in most experimental 
setups. Nevertheless, there is a fast growing literature 
in suspended graphene samples 2 that will eventually tell 
us much more about electron-electron interactions in this 
amazing material. 

In graphene, the strength of the Coulomb interaction 
relative to the kinetic energy is given by the dimensionless 
coupling constant (also called graphene's "fine structure 
constant" ) , 
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where v is the Fermi velocity, and we absorbed the en- 
vironmental dielectric constant (eo) into the definition of 
the charge e. For graphene in vacuum, as in the case of 
suspended samples, a reaches its maximum value, a ~ 
2.2, i.e. interaction effects are expected to be strong. One 
of the ways strong-coupling effects can manifest them- 
selves theoretically is through spontaneous generation of 
a mass m (chiral symmetry breaking). In solid state lan- 
guage, mass generation is equivalent to the opening of a 
gap A = 2\m\ in the electronic spectrum. Hence, in this 



work, the terms "mass" and "gap" are interchangeable. 
In relativistic quantum electrodynamics (QED) in two 
space (plus one time) dimensions, QED2+1, the study of 
this phenomenon started quite a while ago^ and is still 
going strong today. Graphene is actually different from 
QED2+1 because only the fermions are confined to a 2D 
plane, while the field lines extend through the whole 3D 
space. In addition, the Coulomb interaction in graphene 
can be considered instantaneous since the speed of light 
c is much larger than the Fermi velocity (v ~ c/300). 
Hence, Lorenz invariance is not respected, which reflects 
the non-relativistic, purely band origin of the Dirac quasi- 
particles. For the case of graphene, mass generation has 
been predicted 5 -^, although no experimental signatures 
have yet been detected. An in-plane magnetic field also 
favors an excitonic condensate (gap)£. There has been 
a surge of recent numerical activity on the problem of 
mass generation in graphene (without external field), 8 
and a consensus seems to have been reached that above 
a c w 1.1, mass generation occurs^. 

Another way to analyze this phenomenon is as a func- 
tion of N, the number of fermion species, which for 
graphene is N = 4 due to the spin (2) and valley (2) 
degeneracies. In the strong coupling limit a — > 00, gen- 
eration of mass occurs below a critical N which was esti- 
mated to be N c 7 — 9 5,8 . In experiments, a detectable 
gap has so far been observed only in a situation when 
it is actually due to external factors, such as the pres- 
ence of a substrate with specific symmetry, creating sub- 
lattice asymmetry in the graphene plane and thus mak- 
ing the graphene electrons massive (gapped)^. However, 
it is quite possible, as already mentioned, that in "sus- 
pended" graphene, whose exploration has just begun, the 
gap generically exists due to the strong quasiparticle in- 
teraction. 

Whether graphene breaks overall parity (sublattice 
symmetry) in the process of spontaneous mass genera- 
tion depends on the details of the interactions. Long- 
range Coulomb interactions 5 favor equally parity-even 
and parity-odd combinations of masses in the two Dirac 
cones (valleys). On the other hand in QED2+1 it is usu- 
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ally argued that parity-breaking mass generation is not 
possible, but this has to do with the presence of vector 
interactions in the fully relativistic modet^. 

In this work we study the effect of a finite gap on the 
quasiparticle interactions (in particular modifications of 
Coulomb's law), and the renormalization of the quasipar- 
ticle parameters, such as the Fermi velocity and the gap 
itself. Our main goal has been to compute corrections to 
those quantities under the assumption that the system 
is already massive, e.g. due to external factors explicitly 
breaking the sublattice symmetry, as mentioned above. 
However at the end of the paper we also present estimates 
how the new physics we find can possibly affect the spon- 
taneous formation of a gap via the excitonic mechanism. 
We do not address the issue of parity breaking since only 
single valley Coulomb interactions are considered. 

For massless graphene, the large-N limit was studied 
recently in Refs. [Tll[l2 |. extending earlier results^. The 
present work can be viewed as an extension of those stud- 
ies to the massive case. We also find that for massive 
Dirac fermions unconventional modification of the inter- 
action vertex and fermion's properties can occur, not pos- 
sible for strictly massless quasiparticles. In particular, we 
find that there is a crossover regime for Na 3> 1 where 
the 3D Coulomb law, V(r) ~ 1/r, is modified, due to the 
confinement of the electric field lines, to a 2D Coulomb 
law, V(r) ~ ln(l/r), with strong renormalizations of the 
quasiparticle properties. In this non-perturbative regime 
the photon field is confined to the graphene plane lead- 
ing to a situation where the Dirac electrons form a 2D 
"relativistic" Coulomb gas^. Such electronic state has 
never been observed in nature before and, perhaps, with 
developments in the control of the structure of graphene 
samples, it can be studied soon. 

The paper is organized as follows. In the next section, 
for completeness, we analyze the weak coupling regime 
of small a, which can also be relevant to real situations 
(when the substrate screening is strong) . In section IIIII 
we study the l/N expansion for the electronic proper- 
ties of gapped Dirac fermions and show the existence of 
this new intermediate regime where weak confinement of 
the electric field lines leads to strong renormalization of 
electron-electron interactions. We also calculate the im- 
plications of this new regime in the quasiparticle prop- 
erties. Section llVI contains our conclusions, and in Ap- 
pendix A some estimates related to the excitonic gap 
formation are summarized. 



II. INTERACTION POTENTIAL: 
WEAK-COUPLING REGIME 

Our starting point is a model of two-dimensional mas- 
sive Dirac fermions with a gap A = 2\m\. The Hamilto- 
nian of the system iai 

H = Y,^U v(T -P ±m(J ^p + H I , (2) 
p 



where Hi is the quasiparticle interaction 

Hi = -^n p F(p)n_ p , n P = ^* q+p *q' (3) 
p q 

and the potential V(p) will be specified later. We work 
in a two-component representation, so that <7j,i = 1,2,3 
are the Pauli matrices; do = I is the 2x2 identity matrix, 
often omitted for simplicity, and the vector er = (cti, 02). 
Thus the Hamiltonian ([2]) describes the physics in a sin- 
gle Dirac cone (valley). The two valleys in graphene are 
connected by time reversal, which translates into oppo- 
site signs of the mass term: the sign " + " in ^ corre- 
spond to one of the valleys, while the sign " — " to the 
other one. All formulas that follow are invariant under a 
change of the mass sign. The valley (and spin) indexes 
in @ , ([2]) are omitted for simplicity, and it is understood 
that the spinors and the density n p describe a given 
single valley. 

The low energy electronic spectrum (dispersion close 
to the Fermi energy), for Hi = 0, is given by: 

E£ = ±y / v 2 k 2 +m 2 . (4) 

The interaction Hj renormalizes both v and to, as we 
show later. We will only analyze the case when the sys- 
tem is an insulator, i.e. the chemical potential is in the 
gap (e.g. fixed to be zero). It is quite remarkable that in 
recent experiments in gapped samples the chemical po- 
tential can actually be moved from the electron to the 
hole side through the gap, thus causing a metal-insulator 
transition^. In addition to using h = 1 everywhere we 
also put v — 1 in all intermediate formulas and restore it 
only at the end, when necessary. 

The polarization function, II(q, to), for massive Dirac 
fermions in the random phase approximation (RPA) was 
most recently analyzed in detail in Ref. (and can also 
be deduced by appropriate modification of the Lorentz- 
invariant results in QED2+1-). The result is: 

n< q , u >=-^ + — 

where q is the "3-momentum" , 

q = v/|q| 2 -^. (6) 
In RPA the effective potential is given by: 

1 - Vq (0) n( 9 ) |q| 

Firstly, we study the behavior of the potential for weak 
coupling Na -C 1, and in this part of the work we fix N 
to its graphene value N = 4. The first order correction 
to the static potential is: 5V(q) « (Vq 0) ) 2 II(q, a; = 0). 
After transforming back to real space, we represent the 
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correction by the function C(r), and write the total po- 
tential as 



V(r) sa ^(l + «C(r) + 0(a 2 )). 
By using ((5]), we obtain 

/>oo 

C(r) — — 2(mr) / dxJo(2mrx) 



(8) 



which is shown in Fig. [T] (lower panel) . 

At this point we also find it useful to compare 
our results with the well studied case of conventional 
QED 3+ i— , where interaction effects are also governed 
by a dimensionless coupling, the fine structure constant 
oiqed- Eq.([5]) has the same form (with the substitution 
a — > oiQED = 1/137), leading to the so-called Uehling 
potential. The plot of Cqed(t) in this case is shown 
in Fig. Q] (upper panel}^. It should be mentioned that 
(essentially as a consequence of the uncertainty prin- 
ciple) vacuum polarization effects are expected to be 
strong only at distances below the Compton wavelength 
r < Ac = 1/m. For our case we find, by evaluating the 
long distance asymptotic behavior of the integral ^ 19 , 



C(r > m L ) ~ - J — e" 
mr 



(10) 



At short distances we find: C(r <C m x ) = — tt/2, which 
is the well-known massless limit. 

Thus we observe two main differences between massive 
Dirac fermions and QED: (1) The sign of the correction 
is different, i.e. in graphene vacuum polarization weak- 
ens the potential. This can be traced to the fact that the 
charge itself e 2 is not renormalized in graphene^, while 
it diverges logarithmically at small distances in QED (in 
other words the distribution of the vacuum polarization 
charge is very different in the two cases^). (2) The mag- 
nitude of the correction in graphene is more than 10 
times larger (for typical distances), compared to QED 
(see Fig. [1]). This seems to be related to the dimension- 
ality of the problem. Thus, we conclude that (at weak 
coupling) polarization effects can be appreciable in a wide 
range of distances, especially if one assumes that Eq.© 
can be used also away from its strict applicability limit 
ci«l. Typically such potential modification effects are 
important for calculation of localized energy levels (in 
the gap); so far such studies have not been performed 
experimentally in graphene. 

The renormalization of v and m in the weak-coupling 
regime was addressed in Ref. [l6[ and will not be repeated 
here; the main difference from the massless case is that 
the mass provides an effective infrared cutoff where the 
renormalization stops, and the mass itself increases log- 
arithmically (to first order in a). 
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FIG. 1: (Color online.) Plot of C(r) as defined by Eq.flHJ, for 
the case of graphene (blue line) and 3D QED (red line). 



III. 1/N EXPANSION 

Now, we analyze the limit Na ^S> 1, and proceed to 
evaluate the renormalization of the potential and the 
quasiparticle properties. We view the calculation as a 
two-step procedure, which is not exact but will sim- 
plify the problem technically: first we calculate diver- 
gent terms that originate from intermediate integration 
in the high momentum region (m <C |q| <C A), where 
A is the ultraviolet cutoff for graphene (A ~ 1/a, where 
a w 1.42 A is the lattice spacing) and after that, as a 
second step, we concentrate on renormalization from low 
momenta |q| < m (such renormalization will be more 
severe in the strong coupling limit a — > oo). 



A. High momentum regime: m <C |q| < A 

To implement the first (large momentum) part of the 
scheme, one naturally expects the mass term to be unim- 
portant; more formally, we can expand: 



n(q,w) 



4q 



9 >TO. (11) 



Keeping only the first term, one arrives at the effective 
potential, identical to the one found for the massless 



cast 



.11.13. 



2na 

where we use the notation 



7[fl;|q| 

8 q 



(12) 



(13) 



and q is defined by Eq.©. 

Now we evaluate the self-energy correction to one loop, 
in the limit j > 1, where the second term in (|12p domi- 
nates. The self energy is proportional to 1/N in this case. 
The Green's function is 

G _1 (k, u) = lu — v er-k— m 03 — S(k, uj)+ir]sign(Ld) , (14) 
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where the self-energy at one-loop level is 
cPpde 



(27T) 



G (k + p,^ + e)y(p,e). (15) 



Here Go is the free Green's function. By expanding 

E = cjSo + tNT ■ kE„ + maz^m , (16) 

we find 

Z 

^ ' ^ w — Z(l + Ti v )v<t ■ k — Zm(l + E m )cr3 ' 
where Z is the quasiparticle residue: 



Z- 1 = 1-E . 



(18) 



The calculation of the velocity renormalization is then 
practically identical to the massless caseii, except one 
should keep in mind that the mass provides an effective 
infrared cutoff in the integrals (due to the finite value 
of the dispersion at zero momentum ((H)). Thus we put 
the mass to zero to avoid lengthy formulas and keep the 
above in mind. One finally obtains (Sv(k) stands for the 
velocity correction at one-loop order) 



Sv(k) 



Eo + E„ : 



.16 

'N 



pdp 

W) 



(III; 

W) 



2\-3/2 



where k is the external momentum (neglected in the 
Green's function), and written as an effective infrared 
cutoff (thus we assume k 3> to; in the opposite limit m 
is the cutoff). By performing Wick's rotation lj — > iuj 
(which avoids crossing any poles)^, one finds 



fo(fc)/„ = __ 



ln(A/fc), m<fc<A. (20) 



As expected, the result is this region is identical to the 
massless caseii. 

For the mass renormalization we have, written in more 
detail. 



8m(k) 



m 



So 



.16 

''N 



A 



pdp f did yj p 2 — LU 2 



, (2*) 
1 



(2tt) 



V + oj 
— 1 

(p 2 — UJ 2 ) 2 p 2 — LJ 2 



(21) 



Here the first term in the curly brackets corresponds to 
So and the second one to S m . The final result is 

16 1 

6m(k)/m = — ^— ln(A/fe), m<fc<A. (22) 
7H N 

Finally, the quasiparticle residue Z is determined by 
the behavior of Eo. However in the limit j > 1 one en- 
counters some complications, because even the frequency 
integral in the expression for E ( see (EU) is logarithmi- 
cally divergent (although in the sum Eo + E m this ad- 
ditional divergence is not present). This means that we 



have to use the full, finite g potential from ([12")) 
forming the calculation one finds 



Per- 



a f°° 1 - x 2 

E = - ln(A/fc) / dx , 

7T 1 1 'Jo ( 5 7r/8 + VTT^)(l+x 2 )3/2 

(23) 

and therefore, 
8 1 

Z » 1 7 — ln(W4) ln(A/fe), g > 1, m<fc<A. 

7T N 

(24) 

The results for the mass and velocity renormalization 
([20| . I|22[) can be used to form renormalization group 
(RG) equations for these quantities. The reasoning is 
similar to the one presented, for example, in ref. jll| 
for the massless case. One integrates out the high mo- 
mentum degrees of freedom, i.e. momentum regions 
A > |p| > Ai, and the results vary with the quan- 
tity ln(A/Ai) = I. As evident from ([20]) and (HU) the 
renormalization should stop at a scale ~ m. For m large 
enough and Na 3> 1, the functional form of the potential 
(|T2|) is not significantly affected by the RG flow before it 
stops. The RG equations in that case are: 



dv 

m 

dm 
~df ' 

that have the solutions: 



Nir 2 
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ni(k) = i" ( — J " , v(k) = v (— 



(25) 



(26) 



which are valid in the region m <ti k <C A. Here m, v 
are the corresponding quantities at the ultraviolet scale 
A, i.e. their initial band values at the lattice scale. 



B. Low-momentum regime: |q| <g m 

Now we proceed to analyze the low-momentum region. 
In this limit, the polarization ([5]) can be expanded to 
give: 



n(q,w) 



Ajqp 

12nm 



10m 2 



q < to . (27) 



We keep only the first term, as the other terms decrease 
quite fast in powers of q/m. Notice also that in this 
limit n(q, uj) becomes frequency independent. The cor- 
responding RPA effective potential is: 



nq) 



2na 



|q| +g\q\ 2 /m ' 
where we have defined: 

g = g/6 = Na/6 . 



|q| < m , 



(28) 



(29) 
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a) 



b) 



FIG. 2: Two diagrams contributing to the self-energy at two- 
loop level. The wavy line stands for the potential 



By direct numerical evaluation of the polarization 
bubble^, we actually find that the above formula is valid 
even up to |q| ~ to. 

In the strict long-distance limit |q| — > the above po- 
tential tends to the pure Coulomb potential. However, 
in the limit g — > oo there is an intermediate window of 
momenta, m/g <C |q| < to, where the potential crosses 
over to the 2D Coulomb's law, 



12ttto 1 



m/g <C |q| < to. 



In real space we have 
6 



V(r) 



to In 

N \mr 



(-)' 

\mr J 



9 1 
— » r > — 

to to 



(30) 



(31) 



We also comment on some other situations in electro- 
dynamics when the Coulomb potential can be strongly 
modified. It is known that in compact QED 2 +i linear 
confinement can occur due to non-perturbative instanton 
effects 2 ". It is also possible to have a confining potential 
in ferroelectrics where compact field configurations are 
favored leading to linear confinement of charges 2 - . Inter- 
mediate logarithmic behavior (in real space), similar to 
(|3Tj) , can occur in thin films (of thickness d) with large 
dielectric constant k 3> 1—. In that case the logarithmic 
behavior is limited to the region: d <C r <C nd, i.e. d 
plays the role (formally) of the "Compton" wavelength 
1/to in our case, and nd plays the role of Na/m. It 
was argued that such an intermediate regime can actu- 
ally lead to modification of the variable-range-hopping 
law in systems with strong disorder. Finally, a mecha- 
nism similar to the one found in this work, i.e. due to 
the dominance of fluctuations over the bare potential, 
was explored in the context of high temperature super- 
conductivity models based on QED (where it contributes 
to spinon deconfincment-confinement transition) 23 . 

Now we show that the low-momentum region where 
(|28p is valid, also contributes singularly to the self-energy, 
in the strong-coupling limit a 3> 1. Using again the 
expression (ITS"]) with the potential (|28|) , we obtain for the 
velocity renormalization at one-loop level, 



pdp 



dio 



1 



k (2?r) J (2tt) uj 2 - m- 



;V{\v\) . (32) 



Here we have put the fermion energies E p w to. Because 
the potential is static, there is no residue renormalization 



(Z = 1) in this momentum range, and we have finally 

5v^(k) 3 / 1+g \ 3 . , ... 

— = — In w — hitm/k), 

v N \l + gk/mj N 

m/g <C k <C to . (33) 

Performing the same calculation for the mass, one easily 
finds 



8m {1) {k)/m = 8v {1) {k)/v 



(34) 



Let us also examine the two-loop self-energy which is 
given by the two diagrams of Fig. [2l After some rather 
involved calculations we obtain the first, "rainbow" con- 
tribution (Fig.UJa)) 



-(- 

n 2 v: 



'9 9 \ 



in the low-momentum region m/g <^ k -C to. For the 
vertex correction to the self-energy, shown in Fig. [2Kb) , 
we find: 



-— (oj + ivcr ■ k + 3TOCT3) ln(m/fc) . (36) 
N z 



From here we find the corrections to v, to, Z at this order: 



v 

TO 
Z (2) 



17 1 

"Yiv 2 

25 1 

Ttv 2 
1 



ln(m/fe) , 
ln(m/fe) , 



1 - — \n{m/k) 



(37) 



valid for m/g <ti k <C to. 

Observe that a single logarithm appears at two loops, 
meaning that we do not have a conventional renormal- 
ization group situation (piling up of leading logs) in this 
low-momentum region. This behavior is similar to the 
case of a static Coulomb potential in the massless case, 
where a single log appears up to second order of pertur- 
bation theory2i. 

Let m*,v* be the values of these quantities at the low- 
est scale k ~ m/g <C to. Keeping only the dominant one 
loop contribution, we can estimate 

v*/v = m*/m m 1 + j- ln(.g) + 0(l/iV 2 ) . (38) 

From equations (|4"| . (f3"8"]) one can also deduce the correc- 
tion to the dispersion at low momenta: 



\E k \ PS TO* 



{v*fk 2 
2m* 



1 + ^M 3 ) U+ — 



(39) 

for k ~ m/g <C to. 

On the validity of the above expansions, comparing 

Eq. {33]) and (371) we note tnat in tw0 loo P tne V-^ 
expansion breaks down for N w N c = 17/6, when the 
coefficient in front of the log vanishes and the trend in 
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the renormalization of the velocity is reversed towards 
increasing a. In the strong coupling limit, a — > oo, the 
velocity renormalizes to zero, reinforcing a to be large, 
what indicates the possibility of an instability below the 
critical N. Although this analysis is not directly appli- 
cable to graphene, where N = 4, it is similar in spirit to 
the prediction of a metal-insulator transition in massless 
graphene 5 ^, which would take place around the critical 
value a c ~ 1, where usual perturbation theory in the 
Coulomb potential breaks down. 

In this regard it is also important to explore how the 
change in the potential, leading to the modified shape 
([50)) , can affect the spontaneous formation of a mass gap 
via the non-perturbative excitonic mechanism 5-6 . The 
complete self-consistent examination of this problem is 
well beyond the scope of this work^ 5 -; however some es- 
timates are presented in Appendix A. We can conclude 
(see Eq. (|A3[1 ) that the gap increases as g increases, 



Ae -7rA74+(37r/4) ln(g) = A ~ 



S (3 7 r/4) e - 7 rJV/4 



(40) 



although this increase in small in the perturbative 1/N 
regime when \n(g)/N -C 1. On the other hand when 



N 



ln(g) ~ 1 , 



(41) 



the gap enhancement is substantial. However this strong- 
coupling regime in not accessible within the conventional 
large N philosophy, where g = Na/6 is kept fixed while 
N ^> 1, so that the RPA self-consistent scheme is well 
justified. One can also hope that the results represent 
correctly the behavior for N fixed at its physical value 
(N = 4) with g large at strong coupling (a ^> 1), but 
of course in this case diagrams beyond RPA might be 
important and their influence remains unclear. Therefore 
we cannot make a definite conclusion how the system 
behaves under the condition \n(g) ~ 1, although on the 
surface of things this criterion signifies a transition into 
a new low-energy regime which in itself deserves further 
study. 



IV. CONCLUSIONS 

In conclusion, we have studied the behavior of the in- 
teraction potential and the quasiparticlc properties both 
in the weak-coupling g — Na -C 1 and the "strong- 
coupling" g = Na ^> 1 regimes. In the latter case we 
have found an unconventional regime where the potential 
crosses over from the usual 3D Coulomb potential to a 2D 
logarithmic behavior, as if the field lines were confined in 
the plane. This is due to the fact that vacuum fluctua- 
tions dominate over the bare potential, although they do 
so in a limited momentum range. This physics can also 
lead to unusual renormalization of physical quantities at 
distances well beyond the Compton wavelength 1/m (i.e. 
momenta q <C m), up to distances of order g/m ^> 1/m. 
Such effects have not been studied, to the best of our 



knowledge, in conventional QED due to the smallness of 
ocqed- Iii our case both the mass gap and the velocity 
"keep running" (increase) up to the larger distance g/m. 
Since for graphene in vacuum g ss 9, one can hope that 
this physics is observable. On the other hand, the pres- 
ence of various numerical coefficients makes this some- 
what difficult, e.g. in (|3ip g = g/6 = 1.5 is probably too 
small to create a large-enough intermediate energy win- 
dow. Nevertheless, from purely theoretical perspective, 
the massive case exhibits much richer behavior compared 
to the gapless one. 
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APPENDIX A: ENHANCEMENT OF EXCITONIC 
INSTABILITY DUE TO COULOMB'S LAW 
MODIFICATION AT STRONG COUPLING 

The equation for excitonic pairing leading to genera- 
tion of mass m can be derived by determining the self- 
energy in Eq. (fTS")) self-consistently, i.e. substituting 
Go — * G in that equation. Even if the "bare" mass is ini- 
tially zero, a finite (momentum-dependent) mass is then 
non-perturbatively generated, and obeys the equation: 



m( ) = l f d2k m k IM fc ) 
1 2 J (2tt) 2 ^/jfc 2 + m(jfc) 2 ' 



(Al) 



In the gap equation vertex corrections are ignored (in 
the spirit of the large N approach), and the static limit 
of the RPA potential V(|k| is often used 6 . Naturally, if 
on the right-hand side the mass m(k) — m were the non- 
zero mass already present in the Hamiltonian, then the 
above equation would simply be the first perturbative 
correction to the mass, discussed in Section III. 

It is known that the momentum dependence of the 
mass cannot be ignored, and is important down to mo- 
menta k ~ m(0) when the solution levels off. The full 
analysis results in the presence of a critical N, below 
which the solution exists^. However, in order to esti- 
mate the influence of the potential (|30[) on the gap, we 
use a simplified approach which misses the existence of 
a critical coupling, but provides a good qualitative esti- 
mate (and is similar to the initial attacks on the problem 
of chiral symmetry breaking 3 -). This approach amounts 



7 



to ignoring the momentum dependence of the mass, and 
we then obtain: 

f A kdk / 16 \ m f m kdk /12nm\ m 
m " J m 17 {~Nk ) ]e][\ + J m/ - g 17 {~Nk^J ]E±\- 

(A2) 

Here is given by (U]), and in our units v = 1. The 
interaction potential in the above equation is taken in 
the limit g ^> 1, and we have taken into account the 
fact that the potential shape depends on the momentum 
range, as determined by Eqs. (fT^|) . (|30|) . From here we 
find the solution 

TO « Ke -N/i +{ ^/±) ln(S) = ^(W4) e -T«r/4 (A3) 



It is clear that in the regime ln^A/V <C 1 we have an 
exponentially small solution. Only when ln(g) ~ 1, 
a substantial enhancement of the excitonic gap is possi- 
ble, i.e. the exponential dependence disappears, and the 
gap is proportional to the large cutoff scale A. The con- 
dition jf l n (<?) ~ 1 marks, as expected, a transition to 
a non-perturbative regime, where the perturbative 1/N 
corrections previously computed (see Eq. (|38|) ) become 
large, and the 1/./V expansion breaks down. 
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